Coulomb Drag and Spin Coulomb Drag in the presence of Spin-orbit Coupling 
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Employing diagrammatic perturbation theory, we calculate the (charge) Coulomb drag resistivity 
Pd and spin Coulomb drag resistivity pn in the presence of Rashba spin-orbit coupling. Analytical 
expressions for po and pn are derived, and it is found that spin-orbit interaction produces a small 
enhancement to pd and pn in the ballistic regime while pD is unchanged in the diffusive regime. This 
enhancement in the ballistic regime is attributed to the enhancement of the nonlinear susceptibility 
(i.e. current produced through the rectification of the thermal electric potential fluctuations in the 
passive layer) while the lack of enhancement in the diffusive regime is due to the suppression by 
disorder. 

PACS numbers: 73.40.-c, 73. 21. Ac, 71.70.Ej 
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I. INTRODUCTION 

Recent years have seen the emergence of a strong in- 
terest in exploiting the manipulation of the spin de- 
grees of freedom in solid state systems that could lead 
to enhanced performance in electronic devices; this 
field of spin electronics, or simply spintronics 1 , has 
evolved into an exciting subject in condensed matter 
physics. Among the myriad of spintronics proposals, 
gate-controlled Rashba spin-orbit (SO) coupling has at- 
tracted tremendous attention as it offers the possibil- 
ity of all-electrical spin manipulation without the pres- 
ence of a magnetic field. One such experimentally re- 
alizable heterostructure is a double quantum well gate- 
modulated with Rashba spin-orbit coupling. In this pa- 
per, we study theoretically the effect of spin-orbit cou- 
pling on the Coulomb drag properties of 2D bilayer sys- 
tems. As is well known, there exists Coulomb coupling 
between the barrier-separated two-dimensional electron 
gas (2DEG) layers induced by the momentum exchanges 
of the Coulomb-scattered electrons residing in individual 
layers. This phenomenon of Coulomb drag (in the ab- 
sence of spin-orbit coupling), first observed experimen- 
tally by Gramila et alrk between two 2DEG layers, has 
continued to be a subject of thorough investigation 3 ; with 
a theoretical description developed in RefsA^Sii. The 
physical mechanism at work in Coulomb drag can be un- 
derstood as follows: The application of a current through 
one layer (the active layer) causes thermal fluctuations 
in the electron density, which are transferred across the 
barrier to the other layer (the passive layer) due to mo- 
mentum exchanges by electron-electron scattering. In the 
passive layer, the variations in electrical potential asso- 
ciated with the thermal fluctuations in electron density 
then produce a current due to second-order rectification 
effect, quantified by the nonlinear susceptibility which is 
the response function connecting the random potential 
fluctuations and the induced electrical current. The drag 
resistivity p D , which is defined as p D = E passive /j ac tivc 
- the ratio of the electric field strength developed in 
the open-circuited passive layer to the current density in 
the active layer, is a useful experimental measure of the 



Coulomb drag phenomenon. 

Although the subject of Coulomb drag has existed 
since the first experiment more than a decade ago, it 
remains a very vibrant research topic today: recent 
experimental 8 and theoretical^ studies have revealed that 
bilayer hole drag in the low density regime is a factor 
of ~ 10 2 — 10 3 larger than the corresponding electron 
drag. Moreover, other more recent experimentsifl*ii have 
shown that the low-density bilayer hole drag in the pres- 
ence of an in-plane magnetic field (i.e. magnetodrag) has 
a qualitative dependence on the applied magnetic field 
unexpectedly similar to that of the single-layer magne- 
toresistivity, and this surprising observation has been ac- 
counted for theoretically by properly considering the sup- 
pression of screening by carrier spin polarization within 
the framework of the standard Fermi-liquid theory^. In 
light of these recent developments, a very interesting 
question is whether spin-orbit coupling (which is effec- 
tively a momentum-dependent in-plane magnetic field) 
would change the qualitative behaviour of the drag resis- 
tivity. To our knowledge, the effect of spin-orbit coupling 
on Coulomb drag has not yet been investigated, and this 
problem is highly relevant in view of the current growing 
interest in viable semiconductor spintronic systems and 
especially potential device applications employing gate- 
controlled spin-orbit coupling. 

The bilayer Coulomb drag provides direct quantita- 
tive information about carrier interaction effects since 
the drag voltage in the passive layer arises entirely due 
to the carrier-carrier scattering between the layers. The 
study of bilayer Coulomb drag properties in the pres- 
ence of spin-orbit coupling effects therefore introduces a 
method of investigating the interplay between spin-orbit 
coupling and interaction effects. Since both spin-orbit 
coupling and interaction can be tuned by applying exter- 
nal electric and magnetic fields respectively, and also by 
changing the 2D density, the possibility exists for a de- 
tailed understanding of the interplay between spin-orbit 
coupling and many-body interaction effects through bi- 
layer drag studies. This is a main motivation for our 
theoretical work. 

On the other hand, it has also been proposed that 
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there exists an analogous effect, called the spin Coulomb 
dragi^, in which spin-up and spin-down electrons within 
one single device structure (three- or two-dimensional) 
play the roles of the electrons in each individual layers 
in the charge Coulomb drag problem, and now it is the 
Coulomb scattering between the spin-up and spin-down 
carriers that causes the frictional force, damping the rel- 
ative motion of the two spin components. Therefore, un- 
like ordinary charge current, the flow of opposite spin 
carriers (i.e. spin current) is a non-conserved quantity 
even without spin-orbit coupling and tends to be sup- 
pressed because of this drag effect. In the same spirit 
as p£>, one gains a measure of the spin drag by defining 
Pi I = Ef/ji, the spin drag resistivity given by the ratio 
of the gradient of the electrochemical potential of spin-up 
carriers to the current density of the spin-down carriers, 
with the current of the spin-up carriers held zero. This 
effect has been confirmed in a recent experimental paper 
by Weber et ali^ where it was found that the suppression 
of the measured spin diffusion coefficient relative to the 
charge diffusion coefficient could be explained by a cor- 
rection factor 1/(1 + \pti\/p) coming from the theoretical 
spin drag prediction. We note that while the Coulomb 
drag refers to a bilayer system, spin drag refers to a single 
2D layer (or a single 3D system) with the up and down 
components of the spin playing the role of the "active" 
and "passive" carrier components. 



The purpose of the present paper is to investigate the 
effect of spin-orbit interaction on the drag phenomenon 
in semiconductor structures including the usual (charge) 
Coulomb drag and the spin Coulomb drag. We fol- 
low Refs4£ and apply diagrammatic perturbation the- 
ory to calculate the Coulomb drag resistivity pp> and spin 
Coulomb drag resistivity p-^ in the presence of Rashba 
SO coupling. We find that for clean samples the Coulomb 
drag resistivity and the spin drag resistivity are enhanced 
by a correction factor Apo t n/ Pd,h ~ C(7 2 ) m the pres- 
ence of SO coupling (7 is the dimcnsionless spin-orbit 
coupling strength), whereas for dirty samples the spin- 
orbit coupling correction to the Coulomb drag is essen- 
tially completely suppressed by disorder. 



Our paper is organized as follows: we begin in section 
II with the problem of Coulomb drag, where after reca- 
pitulating a set of useful formulas related to Rashba SO 
coupling and Coulomb drag we proceed to evaluate in 
both the ballistic limit (subsection A) and the diffusive 
limit (subsection B) the nonlinear susceptibility, the cen- 
tral quantity in the problem. These results are then used 
in subsection C to calculate the Coulomb drag resistivity 
in the presence of Rashba SO coupling. With the aid 
of the formalism applied to the Coulomb drag, we turn 
in section III to the problem of spin Coulomb drag and 
calculate the spin drag resistivity in the ballistic regime. 
Finally conclusions are presented in section IV. 



II. COULOMB DRAG 

The electronic Hamiltonian in the presence of the 
Rashba SO coupling 



H 



P 

2m 



+ a (a X p) ■ z, 



has the following eigenstates 
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which is labelled by the quantum number p = ±1 called 
the chirality. The corresponding eigenenergy spectra fol- 
low as e p ± — p 2 /2m =F ap. We note that the transforma- 
tion which diagonalizes the Hamiltonian Eq. Q is a local 
transformation, dependent on the momentum p through 
tan^ = p y /p x : 



V2 



oi<t> 



1 

-ie l - 



(3) 



In the presence of Rashba SO coupling, the Fermi energy 
sf = h 2 k 2 F /2m of the system is lowered in comparison 
to the case without SO coupling for the same number of 
electrons and there exists two branches of Fermi wavevec- 
tors Icf± corresponding to the two chiralities. From 
£fc =F OLhk — ef = we can express the two solutions k± in 
terms of hp] on the other hand by particle conservation 
we have n + + n_ = k^/iir + fc 2 _/47r = kp/2it. Solving 
then yields the Fermi energy and Fermi wavevectors in 
the presence of Rashba SO coupling: 



sf = £fq (1 - 27 2 



k F ± = fcFo(±7 + Vl^ 



r 



(4) 
(5) 



where Efo and kpo ar e the Fermi energy and Fermi 
wavevector without SO coupling, here we also define the 
dimensionless spin-orbit coupling strength 7 = a/vpo <C 
1. The density of states v at the Fermi level is readily 
found as 



v± = v\ l± 



7 



(0) 



where v — m/2i:h 2 is the density of states per spin. In 
the following we proceed to evaluate the Coulomb drag 
resistivity of a double-layer 2DEG system with Rashba 
spin-orbit coupling. Using linear response theory, the 
drag resistivity is obtained from the diagrammatic ex- 
pansions in Fig. 1 and is given bya£: 



PD 



V 



167rfc s rCTicr 2 ^ J sinh 2 (hjj/2k B T) 



Ti x (q : uj)T2 X {q,uj)\Ui2(q,uj)\ . (7) 
where a = 2e 2 vD is the Boltzmann conductivity of the 



3 





FIG. 1: Diagrams contributing to the drag resistivity Eq. £J. 
For Coulomb drag (spin Coulomb drag) , the wavy lines repre- 
sent interlayer (intralayer) Coulomb potential; the vertices on 
the left and on the right denote charge current in layer 1 and 
layer 2 (spin-up current and spin-down current) respectively. 



r (<?^) = -J— [ de[n F (e + u>) -n F (e)] Y]tr 
2m J 

{ [Gf(p -q)- G?{p q)) G^ +u (p)J(p)Gf + Jp)} 
+ {q,u^-q,-uj}, (8) 

where Gf ' A (p) = {e - H ± ih/2r)- 1 denotes the 
retarded/advanced Green function, and tr the trace. 
Note that in the expression of the drag resistivity 
Eq. |0, because of the denomenator sinh 2 (Huj/2kBT), 
the dominant contribution of the s integral comes from 
the region ftw < fc^T, and therefore at low temperatures 
the dominant contribution to the drag resistivity comes 
from small values of lu. In the following sections we eval- 
uate Eq. (JHl in respect to two regimes: ballistic regime 
{lot > 1 or ql > 1) and diffusive regime (u)T,ql -C 1), 
where I — v F o T is the mean free path. 






FIG. 2: The three-point vertex representing the nonlinear 
susceptibility V. The vertex on the right side represent cur- 
rent vertex while the two vertices at the top and the bottom 
represent charge density vertices. 



A. Ballistic limit 

In the ensuing discussions we express all Green func- 
tions and currents in the chiral basis, making use of the 
transformation Eq. J2J. To this end we define the trans- 



formation U-\. 
plicitly, 



and U- 



U 



U}Up. 



(9) 



individual layer (D = epQT/m is the diffusion constant), 
the subscripts T','2' signify quantities in layer 1 and 2 re- 
spectively. U\2 is the screened interlayer potential which 
is obtained from solving the corresponding Dyson equa- 
tion in the random phase approximation (RPA), T is the 
nonlinear susceptibility and is given from the three-point 
vertex diagrams in Fig. 2 as 



where 9 — <fi p + q — <f> p = <f> p — 4> p - q is the scattering 
angle from momentum p to p + q (or equivalently from 
momentum p— q to p). It will be useful to also note that 



U 



since J/+._ is a symmetric matrix. 



After a change of variable in e in the second term and 
a transformation into the chiral basis, Eq. JHJ can be 
written as 



T(q,u) 



2iri 



-u 4 



de [n F (e + «) - n F (e)] £tr {/7- [Gf(p - q) - Gf{p - q)] ulG? + Jp)J(p)G? + Jp 

v 

G? +U (P + ~ Gf + JP + 9)] UlG*(p) J(p)Gf (p)} , 



(10) 



where the tilda represents quantities expressed in the chi- 
ral basis. In the ballistic limit, the charge current in the 
chiral basis is given by 



J, X 



aa 



y 



m 

hk x /m — acos0 iasin0 
— iasnxfi hk x /m + acos</> 



(11) 



In the following we will denote the matrix elements of J x 
as Jij for brevity. Here it will be useful to note that the 
off-diagonal elements are related by the hermiticity of J 
as J 
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J. 



21 



-J: 



21- 



Consider the energy integration of the first term in 
Eq. O: 
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J de [n F (e + w) - n F (e)\ tr {[/_ [df(p - q) - G?(p - q)] Ul_G* +ul {p)j (p)G? + Jp)} 

n F {e p _ q+ ) - n F (e p+ ) \ t n _ ^ ^ J n F (e p _q_) - n F (e 



i Im 



22 



(1 



cos60 Im 



LU + £p-q+ - £p+ + i0+ 

n F (e p - q -) - n F (e p -) 



cp-q- 



(1 



(1 



cos60 Im 



u; + e p -q- - e p + + i0+ 
n f (e p _g + ) - n F (e p -) 



+2tt Ji2sin0Im 



p_ + i0+ J ^ cj 

n F (e p _g + ) - n F (e p+ ) _ n F (e p _g_) - n F (e p _) 
i0+ + £p-q- - e p - + iO 4 



a; - 



fc p— q+ 



tp+ 



Cp— q- 

<^(£p+ 



-ep. 



(12) 



We notice that the last term in Eq. (|12|) contains a delta 
function 5 (e p + — e p -) = S (2ap), which forces this term 
to vanish identically upon performing the momentum in- 
tegration in Eq. (|10|) . and hence, rather unexpectedly, 
the off-diagonal terms in the current Eq. (|llfl drop out 
and only the diagonal terms contribute to the nonlin- 
ear susceptibility. The same analysis above follows for 
the second term in Eq. UlUfl . and the expression for the 
nonlinear susceptibility reduces to (with the notation 

J++,-- = ^11,22): 



n F (ep- qil >) - n F (e pll ) n F (e p ^) - n F (e p+qfl 



(13) 



Here we notice that the nonlinear susceptibility T(q, 10) 
is enhanced in the presence of SO coupling. In the bal- 
listic limit, it should be emphasized that the nonlinear 
susceptibility without SO coupling is directly propor- 
tional to the imaginary part of the irreducible polariz- 
ability generally^*^, independent of the low-energy long- 
wavelength assumption made on u and x. Such a general 
statement is not true with SO coupling present, where it 
is readily seen that Eq. l|13f) is not, despite the similar- 
ity, explicitly proportional to the imaginary part of the 
irreducible polarizability in the presence of Rashba SO 
coupling^: 



Im 



LJ + €r 



(p+qp. + ^ 



Up to this point we have not yet made use of any sim- 
plifying assumptions and Eq. (|13fl for the ballistic case 
is general. Now, we evaluate Eq. I|13fl m the limit of 
low-energy u <C 1 and long- wavelength 1 < 1, and we 
neglect the interband terms (/i 7^ //) as these are readily 
shown to be smaller than their intraband counterparts 
(/j = //) by an order of 0(q 2 /k Fa ). At low tempera- 
tures fc_eT <C £ F o, we expand the Fermi functions in e 
in Eq. (|13|) and keep up to O(ou). The evaluation of a 
typical intraband term then goes as 

E Jn(p) (1 + cos*) Im + )-^p + ) | 

^ \ W + 6p_q+ - £p+ + l0+ J 

(14) 



2uk 



FO 



U + X 



T 2 H 



v/1 - 7 I v/I - 7 



where u — uj/v F oq and x = q/2k F Q, v F q = fikpo/m is the 
Fermi velocity without spin-orbit coupling. Evaluating 
the rest of the intraband terms in the same manner, we 
find the nonlinear susceptibility Eq. (|13|) as 



AvDq u 



£fo \J\ — 7 2 
\vDq ( 1 , 
£fo \ 2 



(15) 



1 



2 E E 

m,m'=± p 



(1 + fifi'cosO) 



n F (e 



(16) 
n F{e p +gn) 



e pt j.' e p+qp. 



i0+ 



The reason here is that the current term Jn,22(p + q) — 
Ju,22{p), with SO coupling (c.f. Eq. ifTTl) '). is no longer 
simply proportional to q and so cannot be taken outside 
of the momentum integral. In the long-wavelength low- 
energy limit, however, proportionality is restored as can 
be seen by comparing with the following expression of 
the irreducible polarizability evaluated in the same limit: 



IL(q,u) ~ 2v 1 



r 



(17) 



Therefore, in the ballistic regime, the enhancement of 
the nonlinear susceptibility Eq. Ijl5|l is physically at- 
tributable to the SO-induced enhancement of the intra- 
band electron-electron scattering. 



B. Diffusive limit 

For low temperatures kgT <C £fo, we expand the 
Fermi functions in e in Eq. ijSJ) and keep up to 0(lu). 
Transforming into the chiral basis, T can then be ex- 
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pressed as 



r(*,«) = -££tr (18) 
v 

[U- [GA(p-q)-G*(p-qj\ Ul_ 



Since the spin-orbit coupling strength is assumed to be 
weak 7< 1, we retain terms up to the first nonvanishing 
order of (D(^ 2 ). In the diffusive limit u<l/T,g<l/l, 
we also expand Eq. l(T%|) in terms of these parameters 



up to the first nonvanishing order, e.g. 
function we have 



G eF (p ± q) 



g £f (p)±<i 



■ pi — 



aa- 



for the Green 



Ge F (p) (19) 



For the charge current vertex in the Rashba model, tak- 
ing account of the diffuson pole correction exactly can- 
cels the spin-dependent term in the charge current, leav- 
ing only the bare current term J x {p) = hk x /m. Using 
U\_ _ = Ul _, the nonlinear susceptibility can be writ- 
ten as (dropping the energy and momentum labels of the 
Green function which are understood here as p and Ef)' 



r x (g,^-^£^tr{{Rc{ 
7rz m, 11 I 



G H - G J 



Ul 



iq ■ pirn 



(G R ) 2 - (G A ) 



A\2 



u 



}}g a g r ] 



(20) 



The first term inside the summation is readily shown to 
be identically zero after the trace is taken, in which case 
we simply have 



T x (q,tu) = > cos 

7T 



<j) p ( — ) tr 
[{G R f - (G A f 



u: 



(21) 

)g a g r ). 



We proceed to evaluate the integral in Eq. (|21|l in the 
following. In the transformation U+ ■ ■ ■ f/l, we expand 
and retain terms up to the first order in q/kpa, i-c 
licos^ ~ 1±(1 — g 2 sin</)/2fc 2 ) and sin9 ~ qsm<fr/k. Phys- 
ically, this means taking account only of the intraband 
contribution but neglecting the interband contribution 
as this is smaller than the interband term by an order 
0(q 2 /k F0 ). Eq. (|2"T|) is then evaluated as 



cos2( ; 

p 

tr • 



hmiu+ [(G R f - (G A ) 2 ] U* + }G A G R } 



— - (V+ + U-) 



Ei, P 
cos <b„a — 
p m 

p 



(22) 



tr {im \u + a z \[G R ) 2 - (G A ) 2 ] C/;} G A G R } 



= 2na ¥ 



kF+ kF- 



(23) 



We note that it is crucial, as emphasized in Ref.- , that 
in evaluating the above integrals Eqs. I|22|) - (|23|l the asym- 
metry between electron and hole spectra has to be taken 



into account, for otherwise electron and hole drag would 
compensate each other completely, yielding a null result. 
Using Eqs. ©-© for k± and u±, we find that the SO 
coupling term Eq. (|23ll vanishes and the nonlinear suscep- 
tibility Eq. I|21|) . in the diffusive limit, remains unchanged 
in the presence of Rashba SO coupling: 



2bjq 



ujq /ry 

7T \hJ 



(24) 



In the above calculation we have only taken into account 
the diffuson vertex correction to the charge current ver- 
tex. In the following we include also the diffuson pole 
correction to the charge density vertex, whereby Eq. i(2"T|) 
becomes 

- -v^ cos ^ p (i) tr { lm { u+ (i- aaz ) 



p 
i 



(G R ) 2 - 



Dq 2 + iuj 



1 



Dq 2 — ioj 



(G 



A\2 



u: 



(25) 



Evaluating the integrals, we find the complete expression 
of the nonlinear susceptiblity in the diffusive limit with 
Rashba SO coupling 



r a (g,w) = -4 



ujDq 2 vDq 



{Dq 2 



Ul 2 £fo 



(26) 



which is unmodified in the presence of SO coupling. This 
rather peculiar result for the Rashba SO coupling can be 
explained in light of the following. First, the diffuson 
vertex renormalization restores the charge current to its 
bare value p/m without SO coupling. Second, it will 
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be recalled that, without SO coupling, the nonlinear sus- 
ceptibility Eq. (|18f) in the diffusive limit is commensurate 
with the imaginary part of the irreducible polarizability 
only to the lowest order of ql and uri^. We find that 
the same statement holds true in the presence of SO cou- 
pling. The polarizability in the presence of disorder is 
determined as 



n 




where it is seen that the leading term is unmodified in 
the presence of SO interaction, and the SO coupling cor- 
rection comes in only in the second order (ql) 2 . It follows 
that the polarizability II(q, oj) and therefore the nonlinear 
susceptiblilty T(q, w) remain unchanged to the lowest or- 
der in ql and cut in the presence of SO coupling. It follows 
that, in the drag resistivity Eq. Q, since the dominant 
contribution to the momentum integral comes from small 
q at low temperatures, pd also remains unchanged in the 
presence of SO coupling, i.e. the SO coupling correction 
to the Coulomb drag becomes completely suppressed in 
the presence of disorder. 



C. Drag Resistivity 

By solving the Dyson equation for the the interlayer 
potential in the random phase approximation, U\2 in 



Eq. 10 is given as 4 
U12 (q,u) 



47re 2 sinh(<7d)rii (q, uj) n 2 (q, lo) ' 



(28) 



where d is the interlayer spacing; II (g, oj) is the polariz- 
ability Eq. 1|16|) for the ballistic case or Eq. I|27|) for the 
diffusive case. In the low-energy long-wavelength limit 
U, x <C 1, the expressions II ~ 2v or II ~ 2vDq 2 / (Dq 2 — 
iui) for the ballistic or the diffusive regime remain un- 
changed in the presence of SO coupling, from which it 
follows that the expression of the interlayer potential 
Eq. (|28l) also remains the same within the assumptions 
made. 

The drag resistivity p^ in the diffusive regime and pb in 
the ballistic regime for two identical 2DEG layers without 
spin-orbit coupling are well knowr 



,4,6,7. 



H IT 2 fkuT 



Pd 



Pb 



12 V £fo 



ln- 



1 



h 7T 2 £(3) (k B T 
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2X (qtfgO (kFoiy 



(qT F d) 2 (k F0 d) 2 



(29) 



(30) 



where <7tf = 47re 2 ^ is the Thomas-Fermi wavenumber 
and ksTo = hDq^/d. In the presence of electron- 
electron interaction, the spin-orbit coupling strength a 
will be rcnormalized, becoming a function of the interac- 
tion parameter r s — \^2e 2 m/h 2 kpo. The detailed evalu- 
ation of the many-body renormalization of the SO cou- 
pling is beyond the scope of this paper, nevertheless we 
note that in the statically screened case the renormaliza- 
tion of the SO coupling a amounts simply to the mass 
renormalizationii a(m/m*), which however may not be 
true if dynamical screening is taken into account. Com- 
bining Eqs. J7J, l|15l) and 1261) . we immediately arrive at 
he results that: (1) the drag resistivity is unchanged in 
he diffusive limit pd = Pd', (2) the drag resistivity in the 
ballistic limit is enhanced as 
(27) 



Pd = Pb 



1 + (7T 



(31) 



where we have denoted the renormalized SO coupling 
strength as 7* = j(m/m*) within the static screening 
approximation made. 



III. SPIN COULOMB DRAG 

In the following we perform a calculation of the spin 
Coulomb drag resistivity in the ballistic limit. Analogous 
to Eq. Q the spin Coulomb drag resistivity is given as 



PU 



IQirksTa^ai 



E 



sinli (hw/2k B T) 



IV (q,u)Ti x (q,w) \U(q,u))\ , 



(32) 



where U is the Coulomb potential, a^i is respectively 
the Boltzmann conductivity for spin up and spin down 
carriers. 

First we consider the case without SO coupling. The 
spin-up and spin-down current are given by 



J\x = 
J[x = 



2m 
2m 



1 





1 



(33) 
(34) 



Using these expressions in Eq. i|l(J|) yields 

^f(ept) - n F (e p+q {) 



qr 
m 



E Im {; 



e pT £ p+<jT 



qr 



= ImITf(g,u;) 



(35) 



and similarly for down spin. Putting into Eq. Q32JI gives 
the expression for the spin drag resistivity without SO 
coupling 



PU 



167re 2 /cB Tn 



E 



sinh 2 (hu:/2k B T) 



ImlTrlnJFjC/l 2 , 



(36) 
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consistent with Refsii^iSii^. We note however that our 
result is smaller by a factor of four because in our defini- 
tion of the spin up and spin down currents Eqs. J^-EHl 
we have included a factor of 1/2 for spin-1/2. 

We now consider spin Coulomb drag in the presence 
of Rashba SO coupling. Since in the ballistic limit the 
spin current with Rashba SO coupling is diagonal in spin 
space J sx = (Hk x /2m)a z , the spin up and spin down 
currents are still given by Eqs. (|33fl - l|34[l . Transforming 
them into the chiral basis and substituting the resulting 
expressions into Eq. i|l(J|) for the nonlinear susceptibility, 
it can be shown that, similar to Eq. (|12p. in the ballistic 
limit the off-diagonal terms of the currents J^ x , J[ x do 
not contribute and the nonlinear susceptibility for T^ x 
or rj, x still maintains the same form given by Eq. (jl3|l . 
where now Jn.22 — hk x /Am. Evaluating the momentum 
integral then gives 



vDq 



tFo (1 - 7 2 ) 3/2 
vDq ( 3 9 



(37) 



For the Coulomb potential, we can use in the low tem- 
perature regime the statically screened expression, since 
by virtue of Eq. I|17|) the real part of the polarizability 
is unchanged by SO coupling, the Coulomb potential is 
simply given by 



U(q) 



2ne 2 



<Ztf 



(38) 



Putting Eqs. into the Eq. we find in a 

paramagnetic system n-f = = n/2 the spin drag resis- 
tivity 



pu - tit 



H 7T 
^3 



3( 7 * 

q 



Qtf 



2k 



dq- 

n (q + q T F) 



H 7T 

?3 



1 - 
In ( 1 



•3(7*)' 

2k F0 
9tf 



/ Qtf 



\2kFc 



k B T 

£F0 



(39) 



1 



1 + gTF/2fc F0 



where in the last line in order to evaluate the divergent 
q integral we have introduced an upper cutoff at q = 
2k fq. Our result Eq. (|39|l is in agreement with Refpi& 
in the absence of SO coupling. We observe here that 
the presence of the SO interaction also enhances the spin 
Coulomb drag. 

We comment in passing that, in the diffusive regime, 
correlated impurity scattering has to be taken into 
account since the spin-up and spin-down electrons are 
scattered by the same set of impurities in the sample, 
as emphasized in Rcf. 1 "^. In this paper we have not 



attempted to include the diffusive limit as the effect 
of correlated impurity scattering requires taking into 
account a different set of diagrams^ and is beyond the 
scope of the current paper. However, taking account of 
such correlated impurity scattering effect is expected to 
further enhance the spin Coulomb drag in the presence 
of SO coupling also. 



IV. CONCLUSION 

Before concluding, we point out that our theory indi- 
cates rather small effect of SO coupling on the Coulomb 
and spin drag, at least in the weak SO coupling regime 
considered in our work. The drag is either unaffected by 
the SO coupling or is affected only in the quadratic or- 
der of the SO coupling strength. So the signature of SO 
coupling may not be easy to discern in experimental drag 
measurements. On the other hand, we do find that the 
SO coupling enhances the bilayer Coulomb drag by the 



factor 



i+(7*r 



in the ballistic regime, indicating that 

Coulomb drag would be enhanced by the presence of the 
SO coupling. This is actually consistent with the avail- 
able experimental findings. In particular, the measured 
bilayer hole drag&iS seems to be larger than the corre- 
sponding theoretical results*^^ even after the inclusion 
of many-body effects. Since the SO coupling is strong 
in GaAs holes, the dimensionless SO coupling strength 
(7*) 2 could be large, particularly at the low experimen- 
tal hole densities (because the definition of 7 involves the 
Fermi velocity in the denominator). This could lead to 
an appreciable enhancement of the 2D Coulomb drag for 
low density holes in the ballistic regime bringing experi- 
ment and theory closer together. Unfortunately, we can- 
not make any quantitative comments on this issue since 
our theory is explicitly a weak SO coupling expansion. 

In summary, we have calculated the Coulomb drag re- 
sistivity of a double-layer 2DEG system including spin- 
orbit coupling. In the diffusive limit we find that there 
is no change in the Coulomb drag resistivity while in the 
ballistic limit we find the drag resistivity to be enhanced 

by a factor of 1 + (j*) 2 . We also apply the formal- 
ism to the problem of spin Coulomb drag and find that 
the spin drag resistivity in the ballistic limit is similarly 
enhanced as well. These enhancement effects in the bal- 
listic regime are due to the SO-induced enhancement to 
the nonlinear susceptibility in both the charge Coulomb 
drag and spin Coulomb drag. 
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